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Heat asymmetries in nanoscale conductors: The role of decoherence and inelasticity 

Javier Argiiello-Luengo, David Sanchez, and Rosa Lopez 
Institut de Fisica Interdisciplindria i Sistemes Complexos IFISC (CSIC-UIB), E-07122 Palma de Mallorca, Spain 

We investigate the heat flow between different terminals in an interacting coherent conductor when 
inelastic scattering is present. We illustrate our theory with a two-terminal quantum dot setup. Two 
types of heat asymmetries are investigated: electric asymmetry Ab, which describes deviations of 
the heat current in a given contact when voltages are exchanged, and contact asymmetry Ac, which 
quantifies the difference between the power measured in two distinct electrodes. In the linear regime, 
both asymmetries agree and are proportional to the Seebeck coefficient, the latter following at low 
temperature a Mott-type formula with a dot transmission renormalized by inelasticity. Interestingly, 
in the nonlinear regime of transport we find Ae ^ Ac and this asymmetry departure depends on 
the applied bias configuration. Our results may be important for the recent experiments by Lee et 
al. [Nature (London) 498, 209 (2013)], where these asymmetries were measured. 


I. INTRODUCTION 

Thermoelectrical transport at the nanoscale is a phe¬ 
nomenon of wide interest due to its fundamental and ap¬ 
plied perspectives [I|. From the practical point of view, 
nanostructures spur a wide range of promising thermo¬ 
electric applications such as thermocouples 0, local re¬ 
frigerators Q, thermal transistors [^, and thermal rec¬ 
tifiers among others. The conversion of waste heat 
into electricity seems to be more efficient at the nanoscale 
than at macroscopic scales [^. The fast pursuit toward 
higher efficiency values of the generated electrical power 
in relation to the supplied heat has reached remarkable 
results 0,I1|- However, related fundamental issues such 
as the electronic heat flow traversing a nanodevice still re¬ 
main poorly understood mainly because thermal current 
is not easily accessible in an experiment @ . In many as¬ 
pects, heat flow inherently differs from its electrical coun¬ 
terpart and can reveal information about the number of 
channels available for transport [l^ . the presence of in¬ 
teractions [HI: properties of single-particle wave func¬ 
tions [III , or even superconducting phase differences [l^ ■ 

Recent works have investigated both experimentally 
and theoretically the heat current in atomic-scale junc¬ 
tions Bm. Importantly, power dissipation at atomic 
scales depends strongly on the way in which the transmis¬ 
sion probability varies with energy. Thus, for nanostruc¬ 
tures showing a strongly energy-dependent transmission, 
the measured heat flux is shared quite asymmetrically 
among the contacts whereas for those systems with a 
weakly energy-dependent transmission the heat asymme¬ 
try is strongly suppressed B- This conclusion assumes 
that energy exchange between carriers occurs elastically. 
However, in molecular or atomic junctions, inelastic pro¬ 
cesses can be of critical importance when internal degrees 
of freedom such as rotational or vibrational modes come 
into play B"Bl- As a consequence, these can alter the 
physical scenario. The fundamental question addressed 
in this work is precisely how the heat current asymmetry 
is affected by inelastic processes. We are interested in 
the contact asymmetry Ac, which measures differences 
between the source (J7i) and the drain {J 2 ) heat cur¬ 


rents, and the electric asymmetry Ac, which quantifies 
the heat-current asymmetry in a given electrode when 
the applied voltages Vi and V 2 are exchanged: 

Ac = MVi,V 2) - j2iVi,V2) , ( 1 ) 

Ae = JiiVi,V2) - MV2 ,Vi) . ( 2 ) 

Furthermore, even when only elastic processes are 
present, dephasing mechanisms can also take place. 
Therefore, it is also natural to ask how heat is partitioned 
among the different electronic reservoirs in the presence 
of dephasing. 

To examine both issues, we use the voltage 
and dephasing Bi B probe models, recently general¬ 
ized to treat heat-current flows [271 - B • these formula¬ 
tions, inelastic and dephasing processes are incorporated 
by considering a fictitious terminal attached to the quan¬ 
tum system in such a way that the net electrical and heat 
currents flowing through the probe vanish. In particular, 
for a voltage probe a carrier that enters the probe with 
a given energy is reemitted into the conductor with an 
unrelated energy. In contrast, when only dephasing pro¬ 
cesses are present, the energy-resolved heat and charge 
currents are identically zero at each energy. Since the 
model is independent of the microscopic details of the 
actual scattering mechanisms, the results are simple to 
understand and can be applied to a large variety of sys¬ 
tems. 

Our theory is illustrated with a prototypical model for 
mesoscopic systems: a localized state (representing many 
different quantum systems, i.e., atomic or molecular junc¬ 
tions, quantum dots, etc.) attached to two electronic 
reservoirs and subject to different chemical and temper¬ 
ature biases, as depicted in Fig. [T] 


II. THEORETICAL MODEL 

When a mesoscopic conductor is coupled to i = 1 ■ ■ ■ N 
electronic reservoirs and is driven out of equilibrium by 
electrostatic fields {Vj} or temperature gradients {di}, 
a flow of charge and energy from the reservoirs toward 
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FIG. 1. Schematic of a generic nanoconductor with energy 
level £0 in the presence of a voltage Vi and temperature T<s> 
probes and coupled to left and right reservoirs by tunnel cou¬ 
plings Fi and r 2 . Here, r,j is the tunnel coupling with the 
probe. Vi and T* adjust themselves in order to cancel the 
net flow of heat and charge through the probe. The internal 
potential in the conductor is denoted with U. 


the conductor is established. Charge conservation dic¬ 
tates that all stationary charge flows add up to zero, 
- ~ whereas the sum of thermal currents 

must include the Joule heating term, + 


XiVi) = 0. Within the scattering approach formalism, 
the flows read as [HI 

•Tr, 

f dEA,fj{E), 

(3) 

Ss 

II 

:M 

^ dE {E-fi,)A,f,{E). 

(4) 


The factor 2 originates from spin degeneracy since we 
do not consider external magnetic fields. The heat 
current Ji = — Vili is given by sum of the en¬ 

ergy current jf = i2/h)J2jJ dE{E - Ep)A,^fj{E) 
and the associated Joule dissipating heat power Vjlj. 
The electrochemical potential in reservoir i is defined as 
fii = Ep + eVi with Ep the Fermi energy, and fj{E) = 
[1 -I- exp {{E — is the Fermi-Dirac distribu¬ 

tion function. Each terminal has temperature Ti — 9+9i^ 
obtained from a temperature shift 9i with respect to the 
background temperature 9. The elements Aij = Tr[Jij — 
s\^{E,U(r,{Vk},{9k])sij{E,U(f,{Vk},{9k})] (with k = 

1---N) are given in terms of the scattering matrix s, 
where Tr = Tij is the transmission probability from 
terminal j to contact i and the trace is performed over the 
contact channels. Importantly, due to electronic repul¬ 
sion the potential profile inside the conductor is altered 
when charge is injected by means of electrical or ther¬ 
mal biases. As a result, the scattering properties of the 
conductor expressed by Sij{E,eU{r,{Vk},{9k}) depend 
not only on the carrier energy E but also on the inter¬ 
nal potential landscape W, which depends in turn on the 


set of voltage and temperature shifts. The electrostatic 
response can be determined from the Poisson equation 
e^V^ZJ(r) = —q with £„ the vacuum permittivity and q 
the total charge inside the conductor built up from (bare) 
charges injected by electrical and thermal gradients and 
screened ch arg es created in response to the external per¬ 
turbations [jm.JSl l . 

At sufficiently low biases in the applied voltages and 
temperatures, Eq. ([3]) is expanded up to first order in the 
shifts Vi and 9i and the result can be expressed in matrix 
form: 

(J)^(S^)(I). (^) 

where we have defined the vectors X = [Xi,-- - ,Xn] , 
J = V = [Vip--,VNf and e = 

[ 6 * 1 ,-•• , 6 *Ar]^. The elements of the submatrices G, L, 


M and K are 

the transport coefficients 



JdEiM.S.,-X,){-f'^), 

( 6 ) 

Xy = — 

^ h9. 

l'dE{E-EpmS,,-X,){-f^^), 

(7) 

My = 9Lij 

5 

( 8 ) 

II 

^dX;(A-A^)2(MJy-Ty)(-/'q), 

(9) 


where Aft represents the channel number in the i-th con¬ 
tact and /'q denotes the energy derivative of the Eermi 
distribution function evaluated at Vi = 9i = 0. Equa¬ 
tion ([ 8 |) is a consequence of reciprocity. Additionally, the 
transmission in the linear-response regime is evaluated at 
the equilibrium potential and is thus independent of the 
nonequilibrium screening U. We will later consider the 
nonlinear regime, in which currents do depend on U. 

III. ELASTIC AND INELASTIC PROBES 

To include inelastic processes in the thermoelectric 
transport we consider an additional fictitious probe, de¬ 
noted by d), that plays simultaneously the role of an ideal 
voltmeter and thermometer. Then, both charge X$ and 
heat currents through the probe are identically zero. 
Each current carrier absorbed into the probe is reemitted 
with unrelated phase and energy. We hence use the con¬ 
ditions X$ = = 0 to eliminate the probe voltage Vi 

and temperature and rewrite Eq. with modified 
transport coefficients: 

Gij — Gij J- 

J- — G^^M^j)] , ( 10 ) 

Eij ^ Eij J- E^Gi<^K<^j 

Li<^(^]\At^<^ Et^j )j, (11) 

Kij = Kij J- Z?[Arici>(Af<j,ci>X$j — G$<i>Ar$j) 

J- Miik,{L^^K^j — Ar$$X$j)], (12) 
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and Mij = OLij insofar as the Kelvin-Onsager symmetry 
condition is preserved even in the presence of the probe. 
Here, D = 

When the source of scattering is elastic, one employs 
a dephasing probe. The charge- (heat-) current den¬ 
sity i{E) is determined from the equation 2{E) = 

f dEi{E) [J{E) = f dEj(E)]. We impose the condition 
that for each energy E the probe draws no net current, 
i<s>{E) = j<s>{E) = 0, resulting in the unique probe distri¬ 
bution function fg, = — Substituting /$ 

back into the charge and heat flows, one arrives at 


2e 






A 




Ji = / dE{E - fj,i) ( Aij - 


Aic^A^j 




(13) 

/.■(14) 


Here, Aij — Ai^A,^j /A^^ includes a transmission func¬ 
tion renormalized by decoherence effects due to the probe 
coupling. 


IV. SOURCE-DRAIN CONDUCTORS 


In the following, we focus on a simple geometry: a two- 
terminal conducting device as illustrated in Fig. [T] Let Vi 
(V 2 ) be the bias drop and temperature applied to termi¬ 
nal 1 (2) in the isothermal case (0i = 02 = d)- The power 
measured at each contact is shown to exhibit different 
values depending on the configuration measurement M- 
We commence our analysis with the linear regime in 
which voltage shifts are very small. Due to energy current 
conservation, the condition j7i(Ui,V2) = —J 2 {Vi,V 2 ) 
holds. Hence, the measured heat contact asymmetry 
[Eq. dl])] in the absence of incoherent scattering becomes 

Ac = 2MiiU = -2SGiiVe, (15) 

to leading order in U = Ui — V 2 . Corrections would be 
of the order of U^. In Eq. (fT^ . S = —Ln/Gn rep¬ 
resents the Seebeck coefficient. The asymmetry is pro¬ 
portional to the thermopower 0 since S indeed mea¬ 
sures the asymmetry between electron-like and hole-like 
transport. Importantly, the contact asymmetry [Eq. ([T|)] 
amounts to the electrical asymmetry [Eq. ([2])] due pre¬ 
cisely to the energy conservation condition. Now, in the 
presence of inelasticity (voltage probe) we find that both 
heat asymmetries still coincide (A = Ac = Ac) and are 
given by 

A = 2LiiV9 + 2DV9[G^iL^^Ki^ - G^iK^^Li^ 

— -|- , (16) 


which is valid up to linear order in V. Clearly, when the 
probe is decoupled we recover Eq. m- 

At low temperature, a Sommerfeld expansion of 
Eq. (fTH)) yields 


A = 


3h 



Ti$r$2 


/ 


+ 0{9^), 


(17) 


where the prime indicates that the energy derivative is 
evaluated at E = Ep. Equation (II3 has a surprisingly 
simple form. We recall that in the presence of a voltage 
probe the transmission is split into the coherent term 712 
associated with those carriers that flow between source 
and drain without interacting with the probe and the 
incoherent transmission Ti$T<52/(7i$ -I- 72$), which ac¬ 
counts for the fraction of carriers that are incoherently 
scattered through the probe $ [ 2 ^. Here, we find that 
the heat asymmetry is nicely given by the energy deriva¬ 
tive of both terms summed. In fact, Eq. (HZl can be 
interpreted as a Mott-type formula in which S, which 
is proportional to T [2 in Eq. (fT5|) due to the Mott rela¬ 
tion (39|, becomes modified by the incoherent part but 
keeping the same structural form. 

Eor the dephasing probe, we first perform a linear ex¬ 
pansion for Vij in Eq. 0 . Then, the heat asymmetry 
reads as 

^AJ^jiE(E-Ep (r,, + (-/;,). 

(18) 

An important remark here is in order. The heat asymme¬ 
try for the inelastic probe and the dephasing case differ at 
temperatures higher than the energy scale at which the 
renormalized transmission varies appreciably. However, 
to lowest order in the background temperature, Eq. dUD 
identically gives Eq. (HZD. This implies that at low tem¬ 
perature, the heat asymmetry is largest when the renor¬ 
malized transmission (i.e., the coherent plus the inco¬ 
herent terms) varies rapidly with energy around Ep and 
that both dephasing and inelastic mechanisms contribute 
equally. Deviations appear to higher order in 9. Impor¬ 
tantly, when all transmissions are functions of the local 
density of states [i^, the asymmetry A cancels out in 
the electron-hole symmetry case. 


V. NONLINEAR HEAT ASYMMETRIES: A 
QUANTUM DOT EXAMPLE 

The nonlinear re gim e of thermoelectric transport 
shows unique effects [33, [H, Idll - ld^ . Eor the heat trans¬ 
port, rectifications have attracted a good deal of at¬ 
tention Crucially, electron-electron interactions 

must now be taken into account. It is worthy to men¬ 
tion that a description of the contact heat asymmetry 
in terms of a probe-renormalized transmission is then no 
longer possible. Instead, we need to self-consistently find 
the internal potential of the conductor. For this purpose, 
we illustrate the heat asymmetries for the relevant case 
of a quantum dot coupled to two reservoirs. A quantum 
dot model is the basic description of atomic and molec¬ 
ular junctions in terms of localized atomic or molecular 
orbitals . The scattering matrix is modeled as a Breit- 
Wigner resonance, 




E-eo + iT/2 ’ 


(19) 


ri$ -I- r2$ 
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centered at the atomic/molecular orbital level position 
Eo- Here, denotes the tunneling rates when the lo¬ 
calized level is coupled to the left and right reservoirs 
(Fi = Fi,F 2 ) and Fq = Fi -I-F 2 . The total broadening is 
thus F = Fq -|-F$, where the dot coupling to the fictitious 
probe F$ quantifies the degree of inelastic/dephasing 
processes in our transport description. Due to the sim¬ 
plicity of the Breit-Wigner model, the results for the 
dephasing and voltage/temperature coincide. We leave 
open the question of having different heat asymmetry re¬ 
sponses for more intricate setups. 

The internal potential Li is assumed to be spatially ho¬ 
mogeneous. We thus make the substitution eo —>-£0-1-61^ 
in Eq. (fT^ . U is determined from a discretized ver¬ 
sion of the Poisson equation in terms of a capacitance 
C: CU = qd — 9eq, where qd is the nonequilibrium dot 
charge 



Fi/i+F2/2+F$/$ 
(E - £0 - eW)2 + FV4 


( 20 ) 


and geq follows from Eq. (120)1 by setting all voltages and 
temperature shifts to zero (/i = /2 = /eq and thereby 
/* = /eq)- Note that qd is a nonlinear function of the 
thermoelectric configuration and that U depends implic¬ 
itly on voltage and thermal biases. 

To compute the heat flow in the presence of inter¬ 
actions and inelastic processes, a system of three non¬ 
linear equations are to be solved simultaneously: the 
capacitance equation to obtain L{{{Vk},{9k})) and the 
two conditions for the fictitious probe, 1$ = 0 and 
= 0, that determine Ej, and T$. Then, 
and T,^{{Vi,9i}) are nonlinear functions of the shifts ap¬ 
plied to the electrodes. Once these parameters are self- 
consistently obtained, heat-flow asymmetries can be in¬ 
vestigated. Remarkably and in contrast to the linear 
regime, the contact Ac and electrical A^; asymmetries 
do not generally coincide, 

Ac = 2 J^(Ei, V 2 ) - (El + E2 )I(Ei, E 2 ), (21) 

Ab = J^(Ei,E2)-J^(E2,Ei) 

-EiI(Ei,E2) + E2X(E2,Ei). (22) 


Here, we define J'^(Ei,E 2 ) = Ji^(Ei,E 2 ) = 

—1/2^(El, E2), and X = Ii = —X 2 . Note that Ac depends 
on the particular way in which electrical biases are ap¬ 
plied. For a symmetric electrical bias configuration, Ac 
is indeed a measure of the energy current. Both asym¬ 
metries agree as long as the transmission is symmetric 
under the transformation Ei ^ E 2 , which leads to odd 
charge currents under reversing the bias polarity. How¬ 
ever, this condition is not generally met when interac¬ 
tions are present and rectification effects then arise. Im¬ 
portantly, the Joule heating term affects differently the 
two asymmetries and is, in many cases, the dominant 
contribution, as we demonstrate in the following. 



FIG. 2. Heat-current asymmetry in the linear regime (a) as a 
function of the coupling of the probe F^/Fo, and (b) versus 
the localized level eo/Fo. Parameters: Ep =0, ri = r 2 = 
ro/2, eE = O.OlFo, kgO = O.OIFq. 


VI. NUMERICAL RESULTS 


In this section, we present numerical calculations for 
the heat asymmetries of our two-terminal quantum dot in 
both linear and nonlinear regimes. In either case, the in¬ 
tegrals following from substitution of Eq. (11911 in Eqs. o 
and (HI) require a careful analysis (see Appendix). We be¬ 
gin with the linear regime. In Fig. [2l we show the heat 
asymmetry A = Ac = A^: as a function of the probe 
coupling F$/Fo [Fig.[2Ka)] and the dot level £o/Fo, which 
can be tuned with an external gate voltage [Fig. [2][b)]. 
We observe in Fig. ^a) that inelastic processes reduce 
the heat asymmetry A and that the asymmetry is not a 
monotonic function of the gate. Due to the probe cou¬ 
pling, the dot transmission acquires an additional level 
broadening (we recall that F = Fi-|-F 2 +F$). When 
F$ increases the dot transmission becomes broader and 
shows a weaker energy dependence. As a result, the en¬ 
ergy current is reduced overall. We notice that A shows 
electron-hole symmetry, i.e., A(£o) = —A(—£ 0 ). This 
fact is more explicit in Fig.[5](b) and is due to the absence 
of screening effects in the linear regime. Here, the curves 
A versus the dot level show a resonant-like behavior in 
which A becomes an extremum for 2|£o|/F ~ 1. Addi¬ 
tionally, this A extremal point depends quite strongly on 
kB9 (not shown here). Indeed, at very low temperatures 
the value for which A is maximum or minimum indicates 
the energy scale for which the transmission changes more 
abruptly around the Fermi energy. 

In the nonlinear regime, rectification effects arise, as 
illustrated in Fig. |31 For definiteness, we consider a 


A(2F„ /h) 
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FIG. 3. Heat-current characteristic (^7 — V) for various val¬ 
ues of the probe strength F^/Fo (a) for a symmetric tunnel 
conhguration Fi = 72 = ro/2, and (b) for an asymmetric tun¬ 
nel configuration Fi = O.lFo, and 72 = O.OFq. (c), and (e) 
show enlarged parts oi J' — V corresponding to the quadratic 
and linear behavior of the heat current versus voltage, respec¬ 
tively. (d), and (f) illustrate the contact (Ac), and electrical 
(A_b) heat asymmetries. Parameters: Ef = 0, ksO — O.OIFq, 
eo/Fo = 1, G = 0. 




FIG. 4. Heat-current asymmetry in the nonlinear regime for 
an asymmetric tunnel configuration Fi = O.lFo, and 72 = 
O.OFq. Gontact heat asymmetry Ac = Ji{V) — J 2 {V)-. (a) 
as a function of the coupling of the probe r 4 >/ro, and (b) 
versus the localized level eo/Fo. Electrical heat asymmetry 
A_b = Ji{V) — Ji{—V)\ (c) as a function of the probe coupling 
r.j,/ro, and (d) versus the dot level eq/Vq. Parameters: Ep = 
0, Pi = O.lPo, 72 = O.OFo, eP = Fq, keO = O.OlFo, G = 0. 


symmetrically electrical biased quantum dot, i.e., Vi = 
—V 2 = V/2 with a common background temperature 9 
for both contacts. This electrical and thermal configura¬ 
tion mimics the experimental conditions reported by Lee 
et al. in Ref. 0 - We show the heat flow J = j7i(H) 
through contact 1 for a symmetrically coupled quantum 
dot (Ti = 72 ) in Fig. [31(a), and for an asymmetric tunnel 
configuration in Fig. jSKb) (Fi ^ 72 ). In both cases, we 
observe rectification effects, J{V) yf —J{—V) even for 
moderate voltages. These are mainly caused by the Joule 
heating term, which can be further strengthened by an 
asymmetric potential response in the case Fi 72 [5l|. 
In fact, as shown in Fig. |3Kc), the heat flow becomes a 
quadratic function of voltage, J{V) = 

= Mil represents the leading-order electrothermal 
coefficient M)- Thus, J is quickly dominated by the 
Joule power at low bias Pjouie = oc For a 

small-bias range, Fig. ISTe) displays the linear transport 
regime in which J = oc V (Peltier effect). We 

observe that in the strongly nonlinear regime, the effect 
of increasing F^/Fq [Figs. |3Ka) and |3Kb)] causes a de¬ 
crease of the Peltier and the heat current thus becomes 
more symmetric under reversal of the bias polarity. The 
contact and electric heat asymmetries, Ac and Ac, are 
shown in Figs. ISKd) and (f). We observe that inelas¬ 
tic processes (increasing F^/Fq) reduce the value of Ac 
[Fig. EKd)] since the contact heat asymmetry Ac coin¬ 
cides with the energy current for symmetric biases, 


as shown in Eq. m- Hence, by increasing F^/Fq the 
transmission acquires a weaker energy dependence, lead¬ 
ing to a suppression of the energy current for our device 
and therefore a decrease of Ac . The electrical heat asym¬ 
metry is, by construction, insensitive to rectification ef¬ 
fects, as depicted in Fig.|3l[f). Moreover, we also observe 
a decrease of Ac as the amount of incoherent scattering, 
r<i>/ro, increases. 

Finally, we discuss the behavior of the heat asymme¬ 
tries with the probe coupling strength r$ and dot level 
eo in Fig. [4] for the nonlinear regime (we set eV/Fo = 1). 
The heat-contact asymmetry dependence with F^/Fq is 
presented in Fig. SKa) for different eo/ro values. We 
observe departures from the electron-hole symmetry, 
Ac(eo) 9^ — Ac(—Eo), when the strength of the probe 
is relatively small, whereas if F^/Fq becomes larger such 
effects are removed due to an overall reduction of this 
heat asymmetry. Figure SKc) shows the electrical asym¬ 
metry, which is electron-hole symmetric by construction. 
As previously, Ac is broadly reduced when F^/Fq grows. 
The dot gate dependence of the heat asymmetries, for 
specific values of F^/Fq, is shown in Figs. SKb) and (d). 
In both cases, the heat asymmetries show a peak struc¬ 
ture, which is reduced with increasing probe strengths. 
However, Fig.|4l[b) clearly shows the absence of electron- 
hole symmetry for Ac. 


AGFGh) ^ SA(2r„/h) 
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VII. CONCLUSIONS 


In closing, we have formulated a generic framework for 
the assessment of inelastic and dephasing processes in 
the power asymmetry of nanoscale junctions. We have 
found that in linear response, the heat-current asymme¬ 
tries (both measured in a given contact or in different 
electrodes) agree and are given at low temperatures by 
the energy derivative of a modified transmission function. 
In the nonlinear regime of transport, both asymmetries 
differ and present an interesting behavior in terms of the 
coupling to the dephasing probe and the gate-tunable en¬ 
ergy level. Quite generally, the heat asymmetries vanish 
with an increasing amount of inelasticity or dephasing. 

Our results are independent of the microscopic origin 
of incoherent scattering. Qualitatively, we believe that 
our main conclusions will be robust and applicable to a 
large variety of systems. Yet, it would be highly desirable 
to investigate in future works specific models taking into 
account, e.g., electron-phonon interactions. 

Further extensions of the model should consider cool¬ 
ing effects [^ . whose efficiency in the nonlinear regime 
and in the presence of incoherent scattering remains an 
open issue. Another interesting question, perhaps more 
fundamental, is the development of magnetic-field asym¬ 
metries in multiterminal setups [blj. It is well known 
that in the nonlinear regime departures of the Onsager 
reciprocity are quite general (52j |. The role of inelastic¬ 
ity and decoherence is less clear. Finally, we would like 
to mention the exciting possibility of implementing rec¬ 
tifying nanojunctions for energy harvesting [b^. A deep 
study of the combined effect of nonlinearities and inco¬ 
herent scattering would bring the goal of waste-heat-to- 
electricity nanoconverters closer to reality. 
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Appendix: Charge and heat current integrals 


In our numerical analysis, it is worth to calculate the 
integral for the charge current through the source elec¬ 
trode. 




Im 2 : 



FIG. 5. Integration contour in the complex plane. 


and the corresponding heat flux 
= r ! dE{E-^il) 

d Jm. 


FiF 


ir 2 


{E - ey + r2/4 


X [ME)-ME)]. 


(A.2) 


Here, e can be Eq for the linear response or Eq + eU{V, 9) 
in the nonlinear regime of transport, with U evaluated 
self-consistently in terms of the applied voltage V and 
temperature difference 9. 

An analytical solution of Eqs. (lA.lIl and (IA.2I1 can 
be obtained by noticing that the Fermi functions can 
be expressed in terms of the digamma function ^'(z) = 
T'{z)/T{z): 


fyE) = l 
_ 1 
“ 2 


1 — tanh 


( 

V 2/cbT, j 



TT 



(A.3) 


where Wj{E) = [E — M/i'^^sTj). The first (second) 
has singularities at Wj = in (n + i) [wj = —in (n + 4)] 
with n S W. 

Consider now the integral 


1= [ dE T(E)f, (E) = + Is+Ic, 

Jr 


where t{E) 


{E - ef + (r/2)" 


and 


(A.4) 


Ia = t;[ dET{E ), 

^ JM. 


E- 




2nkBTj J 


2nkBTj 


(A.5) 

(A.6) 

(A.7) 
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We compute these integrals using the residue theorem. 
For IA and Iq we choose the upper semidisk of ra¬ 
dius R while for Ib it is convenient to integrate over the 
lower semidisk (see Fig. [S])- In the limit of infinite 
radius {R oo), the integrals along the external paths 
7 ^ vanish and we find 



.€ — fij — zr/2\ 
2TTkBTj ) 

.e — fij + *r/2\ 

2TTkBTj J 


(A. 8 ) 

(A.9) 

(A.IO) 


Here, (z) = 'k ± paths S^, 7 ^ 

are crossed anticlockwise. Js can be obtained analo¬ 
gously to the charge current case: 


Js = — (e - /xi) Im (z) - (z)] 

- 27ri Re (z) - ( 0 )] . (A.15) 


To compute = J^- + J^+, we use the polar represen¬ 
tation z = 


Substituting in Eq. (IA.4I) and defining = 

(3 + mfei ± '&r-) ■ CD becomes 

Ii = -|f^Im[>b{7)-*(7)] . (A.ll) 

where we have used the property 'k (z*) = 4^* (z). 

The expression for the heat current [Eq. (IA.2|) ] is to 
be treated with caution because in this case there arises 
a nonzero contribution from the integration along 7 ^. 
Consider 




'0 (Ae*®Te)^ + r2/4 


X 



^ -Re-±,A _/l -Re-±,. \ 

27rfcBTi J V 2 27rfcsT2 J 


(A.16) 

Let g±{R,0) be the function under the integral of 
Eq. (IA.16F For |z| —>• 00 we use the asymptotic value 
'k (z) ^ ln(z). As a consequence, for every 6 G ( 0 , 7 r) 
one has (;± (A, d) —> iln(T 2 /Ti). Thus, 

R^(x> 


J= [ dz(z-/ri)r(z) [/i{z) -/ 2 {z)] 

JR 

= lim - Js), (A. 12 ) 

it^oo ZTT 


lim J± =7rzln—(A.17) 

i?—boo 


with 


Js= [ dz{z- pi)t{z) [T+ (z) - (z)] 


Collectinff Eos. (lA.lbl) and (|A.17p in Ea. (|A.2(). we finallv 
obtain the analytical expression for the heat current: 

-F 1 dz{z - ni)T{z)['i>f {z)-'Sf (z)] , 

(A.13) 

Ji= (e Mi)Im[4'(2+) 4'(2^)] 

JiR = [ dz{z- Pi)t{z) [T+ (z) - (z)] 

“'t'r 


+ ^Re[vk(z+)-4.(z+)]-^ln|. 

(A.18) 

-k / dz{z- m) t{z) [TC (z) - (z)] . 

“'T'J 

(A.14) 
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